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Abstract: It is shown by both analytical methods and numerical simulations that ex- 
tremely long living spherically symmetric oscillons appear in virtually any real scalar field 
theory coupled to a massless dilaton (DS theories). In fact such "dilatonic" oscillons are 
already present in the simplest non-trivial DS theory - a free massive scalar field coupled 
to the dilaton. It is shown that in analogy to the previously considered cases with a single 
nonlinear scalar field, in DS theories there are also time periodic quasibreathers (QB) as- 
sociated to small amplitude oscillons. Exploiting the QB picture the radiation law of the 
small amplitude dilatonic oscillons is determined analytically. 
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1. Introduction 

Long-living, spatially localized classical solutions in field theories containing scalar fields 



exhibiting nearly periodic oscillations in time - oscillons - |I[-|lf:] have attracted consid- 
erable interest in the last few years. Oscillons closely resemble "true" breathers of the 
one-dimensional (D = 1) sine-Gordon (SG) theory, which are time periodic and are expo- 
nentially localized in space, but unlike true breathers they are continuously losing energy by 
radiating slowly. On the other hand oscillons exist for different scalar potentials in various 
spatial dimensions, in particular for D = 1,2,3. Just like a breather, an oscillon possesses 
a spatially well localized "core", but it also has a "radiative" region outside of the core. 
Oscillons appear from rather generic initial data in the course of time evolution in an im- 
pressive number of physically relevant theories including the bosonic sector of the standard 
model [|r7|]-pO|. Moreover they form in physical processes making them of considerable 
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importance [£^]-[28|. In a series of papers, |12|, [^]-|^, it has been shown that oscillons 
can be well described by a special class of exactly time-periodic " quasibreathers" (QB). 
QBs also possess a well localized core in space (just like true breathers) but in addition 
they have a standing wave tail whose amplitude is minimized. At this point it is important 
to emphasize that there are (infinitely) many time periodic solutions characterized by an 
asymptotically standing wave part. In order to select one solution, we impose the condition 
that the standing wave amplitude be minimal. This is a physically motivated condition, 
which heuristically should single out "the" solution approximating a true breather as well 
as possible, for which this amplitude is identically zero. The amplitude of the standing 
wave tail of a QB is closely related to that of the oscillon radiation, therefore its computa- 
tion is of prime interest. It is a rather non-trivial problem to compute this amplitude even 



in one spatial dimensional scalar theories |32], |3^. In the limit when the core amplitude 
is small, we have developed a method to compute the leading part of the exponentially 
suppressed tail amplitude for a general class of theories in various dimensions |^l|] . 

In this paper we show that oscillons also appear in rather general (real) scalar field 
theories coupled to a (massless) dilaton field (DS theory). Dilaton fields appear naturally 
in low energy effective field theories derived from superstring models |33, 34, |3^ and the 
study of their effects is of major interest. As the present study shows, the coupling of a 
dilaton even to a free massive scalar field, referred to as the dilaton-Klein-Gordon (DKG) 
theory, which is the conceivably simplest non-trivial DS theory, has some rather remarkable 
consequences. This simple DKG theory already admits QBs and as our numerical investi- 
gations show from generic initial data small amplitude oscillons evolve. We concentrate on 
solutions with the simplest spatial geometry - spherical symmetry. We do not think that 
considering spherically symmetric configurations is a major restriction since non-symmetric 
configurations are expected to contain more energy and to evolve into symmetric ones p5| . 
The dilatonic oscillons are very robust and once formed from the initial data they do not 
even seem to radiate their energy, hence their lifetime is extremely long (not even detectable 
by our numerical methods). 

Our means for constructing dilatonic oscillons will be the small amplitude expansion, 
in which the small parameter, e, determines the difference of oscillation frequency from 
the mass threshold. The small amplitude oscillons of the DKG theory appear to be stable 
in dimensions D = 3,4, unstable in D = 5, 6, and their core amplitude is proportional to 
e^. This is to be contrasted to self- interacting scalar theories whose oscillons are stable in 
D = 1,2, unstable in D = 3, and their core amplitude is proportional to e. The master 
equations determining oscillons to leading order in the small amplitude expansion turn out 
to be the Schrodinger-Newton (SN) equations. The main analytical result of this paper is 
the analytic computation of the amplitude of the standing wave tail of the dilatonic QBs 
for any dimension D, and thereby the determination of the radiation law and the lifetime 
of small amplitude oscillons in DS theories. The used methods have been developed in 



Refs. [H, H, 13, m and 

The above results, namely the stability properties and the SN equations playing the 
role of master equation, show striking similarity to those obtained in the Einstein-Klein- 
Gordon (EKG) theory, i.e. for a free massive scalar field coupled to Einstein's gravity. 
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where also stable, long living oscillons (known under the name of oscillating soliton stars, 
or more recently oscillatons) have been found and investigated in many papers [38|-|44]. 



2. The scalar-dilaton system 



The action of a scalar-dilaton system is 
A = I dtd^x 



-2Kip 



urn 



(2.1) 



where if is the dilaton field and $ is a scalar field with self interaction potential U{^). 
The energy corresponding to the action (|2.1| ) can be written as 

1 



E 



d^x£, 



£ 



2 L 



+ e 



f/($). 



2.2 



where £ denotes the energy density. In the case of spherical symmetry 



E 



dr 



vr 



D/2 D-l 



'0 r(z)/2) 

We assume that the potential can be expanded around its minimum at $ 



(2.3) 



as 



k + 1 



(2.4) 



k=l k=l 

where are real constants. For a free massive scalar field with mass m the only nonzero 
coefficient is gi = m?. If g2k = for integer k the potential is symmetric around its 
minimum. In that case, as we will see, for periodic configurations the Fourier expansion of 
$ in t will contain only odd, while the expansion of (p only even Fourier components. For 
spherically symmetric systems the field equations are 



D 

+ -TT^ + — 



aV D 



r or 



1 dip 
dr 



-2k e 



[/(cK). 



(2.5) 
(2.6) 



Qf2 Qj.2 J. 

Since gi = m is intended to be the mass of small excitations of $ at large distances, we 
look for solutions satisfying — > for r — > 00. Finiteness of energy also requires $ — > 
as r — > 00. Rescaling the coordinates as t — > t/m and r — > r/m we first set gi = vn? = 1. 
Then redefining (p — > ip/{2K,) and ^ ^/{2k) and appropriately changing the constants 
gk we arrange that 2k = 1. If for some reason we obtain a solution for which ip tends to a 
nonzero constant at infinity then the dilatation symmetry of the system allows us to shift 
93 and rescale the coordinates so that it is transformed to a solution satisfying p ^ for 
r — > 00. 

An important feature of a localized dilatonic configuration is its dilaton charge, Q. It 
can be defined for almost time-periodic spherically symmetric configurations like oscillons 
as: 

99 « Q r^--^ for r ^ 00 in / 2 (2.7) 
</? « Q Inr for r — > 00 in Z) = 2 . (2-8) 
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3. The small amplitude expansion 



In this section we will construct a finite-energy family of localized small amplitude solutions 
of the spherically symmetric field equations (2.5) and ( |2.6[ ) which oscillate below the mass 



threshold Ip™. It will be shown that such solutions exist for 2 < D < 6. The subtleties 



of the case D = 6 will be dealt with in subsection 3.5. The result of the small amplitude 
expansion is an asymptotic series representation of the core region of a quasibreather or 
oscillon, but misses a standing or outgoing wave tail whose amplitude is exponentially small 
with respect to the core. The amplitude of the tail will be determined in section ^ 

We are looking for small amplitude solutions, therefore we expand the scalar fields, (p 
and <I>, in terms of a parameter e as 

oo oo 

ip = Y,^'^k, $ = 5]e'=$fc, (3.1) 

k=l k=l 

and search for functions (pk and tending to zero at r — > oo. The size of smooth 
configurations is expected to increase for decreasing values of e, therefore it is natural to 
introduce a new radial coordinate by the following rescaling 

p = er. (3.2) 

In order to allow for the e dependence of the time-scale of the configurations a new time 
coordinate is introduced as 

T = uj{e)t. (3.3) 

Numerical experience shows that the smaller the oscillon amplitude is the closer its fre- 
quency becomes to the threshold uo = \. The function ijo{e) is assumed to be analytic near 
a; = 1, and it is expanded as 

oo 

u:\e) = l + Y,e^i^k. (3.4) 

k=l 

We note that there is a considerable freedom in choosing different parametrisations of 
the small amplitude states, changing the actual form of the function io{e). The physical 
parameter is not e but the frequency of the periodic states that will be given by w. After 



the rescalings Eqs. (p. 5]) and (2.6) take the following form 



, , 2^-1 dip 




a# + ' -9^ + ' "—^ -^^= I + >^ " 1 • (3.6) 



Substituting the small amplitude expansion (3T) into (3^) and (|3.6D , to leading e 
order we obtain 
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Since we are looking for solutions which remain bounded in time and since we are free to 
shift the origin r = of the time coordinate, the solution of ( |3.7D can be written as 



^iiT,p) = Pi{p)cosT , ipi{T, p) = pi{p) , (3.8) 

where Pi (p) and pi (p) are some functions of the rescaled radial coordinate p. 
The terms in the expansion of (|3.6D yield 



^ = [1 + cos(2r)] . (3.9) 

This equation can have a solution for ip2 which remains bounded in time only if the time 
independent term in the right hand side vanishes, implying Pi = and consequently 
$1 = 0. Then the solution of (|3.9D is ip2{T,p) = P2{p)- The terms in (3.5) yield 



^2$ 



2 +$2 = 0. (3.10) 



Since = we are again free to shift the time coordinate, and the solution is <I>2(t, p) 

P2{p) COST. 

The order terms in the expansion of (p.6|) give 



■9^3 _ d'^Pi , D - 1 dpi 

dr^ ~ dp^ ^ p dp- ^'^ ' 

In order to have a solution for Lp^[T,p) that remains bounded in time, the right hand side 
must be zero, yielding pi{p) = pu +Pi2p'^^^ when D ^ 2 and pi{p) = pu +pi2^np for 
D = 2, with some constants pn and pi2. Since we look for bounded regular solutions 
tending to zero at p — > oo, we must have pn = pi2 = 0. As we have already seen that 
$1 = 0, this means that the small amplitude expansion ( |3.1D starts with e"^ terms. The 
solution of ( p. 11 ) is then ip^{T,p) = Psip)- The order terms in the expansion of (^ 



give 



+ $3 - W1P2COST = 0. (3.12) 



9r2 

This equation can have a solution for <I>3 which remains bounded in time only if the res- 
onance term proportional to cost vanishes, implying u>i = 0. After applying an order 
small shift in the time coordinate, the solution of (|3.12| ) is <I'3(r, p) = Psip) cost. Contin- 
uing to higher orders, the basic frequency siuT term can always be absorbed by a small 
shift in T. It is important to note that after transforming out the siuT terms no sin{kT) 
terms will appear in the expansion, implying the time reflection symmetry of $ and ip at 
T = 0. 

3.1 The Schrodinger-Newton equations 

The terms in the expansion of (|3.5| ) and (^]^) yield the differential equations 



5t2 



+ ^>4 



d^P2 ^ D-l dP2 , , ^ 



COST - ^52^1 [1 + cos(2t)], (3.13) 

«9V4 d'^p2 D-ldp2 1 2i, , M f1tA\ 

-TT^ = —nr -, — ^ T "2 1 + cos(2t) . (3.14) 

dr^ dp^ p dp 4 ^ ^ ^ ^' ^ ^ 
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The function $4(r, /o) and ip4{T,p) can remain bounded only if the cost resonance terms 
in (3.13) and the time independent terms in ( p.l4 ) vanish, 



d'^P2 D-1 dP2 



P dp 



dp'^ 



p dp 4 



(3.15) 
(3.16) 



Then the time dependence of <I>4(r, p) and ipi{T,p) is determined by ( 3.13 ) and ( 3.14| ) as 



1 1 

$4(t,p) = P4(p)cosr + -ff2^|[cos(2r) -3] , Mr,p) = Piip) - —P2{pf cos{2t) . 

lb 

(3.17) 

Here we see the first contribution of a nontrivial U{^) potential, the term proportional to 
§2 in $4. If (and only if) the potential is non-symmetric around its minimum, even Fourier 
components appear in the expansion of 
Introducing the new variables 



S = P2 + W2 



(3.18) 



(|3l5|) and (|1|) can be written into the form which is called the time-independent 
Schrodinger-Newton (or Newton-Schrodinger) equations in the literature: 



d^S D-1 dS 
p dp 



dp^ 
d^ 
dp' 



+ 



D-1 ds 



p dp 



0, 



0. 



(3.19) 
(3.20) 



We look for localized solutions of these equations, in order to determine the core part of 
small amplitude oscillons to a leading order approximation in e. The main features of the 
solutions depend on the number of spatial dimensions D. For D > 6 positive monotonically 
decreasing solutions necessarily satisfy s = S, they tend to zero, furthermore, the Lane- 
Emden equation holds 

d'^s D-1 ds 



dp"' 



+ 



p dp 



0. 



(3.21) 



For D > 6 solutions are decreasing as 1/p'^ for large p, consequently they have infinite 
energy. It can also be shown that solutions of the original Schrodinger-Newton system 
with s S, and a necessarily oscillating scalar field, have infinite energy, hence there is 
no finite energy solution for D > 6. For D = 6 the explicit form of the asymptotically 
decaying solutions of ( 3.21 ) are known 



±5 



24a^ 



(l + aV^) 



,2^2 ' 



(3.22) 



where a is any constant. Since the replacement of $ with — $ and a simultaneous reflection 
of the potential around its minimum is a symmetry of the system, we choose the positive 



sign for S in (3.22). For D = 6 the total energy remains finite. 



-6- 



If D < 6, then localized solutions have the property that for large values of p the 
function S tends to zero exponentially, while s behaves as s ~ sq + sip^~^ for D ^ 2 and 
as s « So + si In/) for D = 2, where sq and si are some constants. Since we are interested 
in localized solutions we assume 2 < D < 6. From ( |3.19| ) it is apparent that exponentially 
localized solutions for S can only exist if s tends to a negative constant, i.e. sq < 0. In this 
case the localized solutions of the Schrodinger-Newton (SN) equations ( |3.19D and ( 3.2C| ) 



can be parametrized by the number of nodes of S. The physically important ones are the 
nodeless solutions satisfying S > 0, since the others correspond to higher energy and less 
stable oscillons. 

Motivated by the asymptotic behaviour of s, if -D / 2 it is useful to introduce the 
variables 

"^l^l- - — A".- (3.23) 

In 2 < D < 6 dimensions these variables tend exponentially to the earlier introduced 
constants 

lim /i = si , lim v = sq . (3.24) 

p—^oo p—^oo 

Then the SN equations can be written into the equivalent form 

1 + 1^5^ ^0. (3.25) 
1 + ^5-0. (3.20) 

(fS D-ldS , ^ n , . 

d^^^^d^^^""^^^^^ (3.27) 

The SN equations ( 3.19| ) and ( p. 201 ) have the scaling invariance 

{S{p), s{p)) ^ iX^S{Xp), \h{\p)) . (3.28) 

If 2 < D < 6 we use this freedom to make the nodeless solution unique by setting sq = — 1. 
At the same time we change the e parametrization by requiring 

uj2 = -l for 2<D <6, (3.29) 

ensuring that the limiting value of if vanishes to order. Going to higher orders, it can 
be shown that one can always make the choice uji = for i > 3, thereby fixing the e 
parametrization, and setting 



to 



\J\-£^ for 2 < L> < 6 . (3.30) 



If D = 6, since both s and S tend to zero at infinity, we have no method yet to fix the 
value of OL in (|]2|). Moreover, in order to ensure that tends to zero at infinity we have 



to set 

^2 = for L> = 6. (3.31) 
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3.2 Absence of odd e powers in the expansion 

Calculating the order equations from (|3.5|) and (3.6) and requiring the boundedness of 



$5 and 995 we obtain a pair of equations for P3 and p-y. 

+ {P2 + UJ2)P3 + P2iP3 + ^3) 



(fPs D-1 dP3 



+ 



P dp 



d^P3 D 
dp^ 



-^dp3 I 
p dp I 



These equations are solved by constant multiples of 

dP2 



P3 = 2P2 + p- 



dp 



P3 + = 2{p2 +U}2) + P 



dp2 

dp 



(3.32) 
(3.33) 

(3.34) 



corresponding to the scaling invariance (|3.28| ) of the SN equations. In D > 2 dimensions 
P3 given by (|3|) tends to CJ3 — 2 for large p. Since we are looking for solutions for which 
ip tends to zero asymptotically, after choosing ^3 = we can only use the trivial solution 
P3 = P3 = 0. The important consequence is that ^3 = ip^ = 0. Going to higher orders 
in the e expansion, at odd orders we get the same form of equations as ( 3.32| ) and ( 3.33| ), 
consequently, all odd coefficients of and (fk can be made to vanish. Instead of the more 
general form {3A) we can write the small amplitude expansion as 

00 00 
^e^''^2k, ^ = ^' 



(3.35) 



k=l 



k=l 



3.3 Higher orders in the e expansion 



The order equations, after requiring the boundedness of $6 and tpQ, yield a pair of 
equations for P4 and P4: 



d^P 



dp 



1 h {P2 + i^2)Pi + P2{Pi + i^i) 

p dp 



d^P4 D 
'd^^~ 



P 



1 

~ 32 
1 dpi 
dp 



pi + 



5 2 

792 



793 



Pi 



0. 



+ -P2P4 - 7P2P2 



0. 



(3.36) 
(3.37) 



This is an inhomogeneous linear system of differential equations with nonlinear, asymp- 
totically decaying source terms given by the solutions of the SN equations. Since the 
homogeneous terms have the same structure as in (|]3|) and (Q), one can always add 
multiples of 

' - 2P2 + p—r- , p\' +uji = 2{p2 + UJ2) + p — 



p 



, ) i-^A I ^1 — "vi^z I I H 7 ) (3.38) 
dp dp 

to a particular solution of ( 3.36| ) and ( 3.37 ). If 2 < D < 6 we are interested in solutions for 
which at large radii P4 decays exponentially, and P4 ~ go + Qip"^'^ with some constants qq 
and qi. We use the homogeneous solution (|3.38 ) to make qo = 0. Since similar choice can 
be made at higher e orders, this will ensure that the limit of 93 will remain zero at /? — > 00. 
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We note that, in general, it is not possible to make qi also vanish, implying a nontrivial e 
dependence of the dilaton charge Q. 

The resulting expressions for the original $ and ip functions are 



$ = e^p^ cos r + <j P4 cos r + -g2P2 [cos(2r) - 3] + e*' <( Pe cos r 



p3 / 

^256 ( ^ ^ y ^2 + 853 ) cos(3r) - 52 



+ 



92 



3P2P4 - {P2 + UJ2)Pi 



dp 



P2Pi-{p2+UJ2)P2 + {^ 

cos(2t) i +0(e^). 



1 



p2 

Pi - -|-cos(2t) 



+^52^1 [9 COST - cos(3r)] \ + 0{e^) , 



dP 

^P2Pa-{P2+UJ2)P2 - { -j^ 



(3.39) 



cos(2t) 
(3.40) 



where the functions P2 and p2 are determined by the SN equations ( |3.15| ) and (3.16), P4 
and p4 can be obtained from ( p. 36 ) and ( 3.37] ), furthermore, the equations for Pg and pe 
can be calculated from the order terms as 



d^P, 



dp- 



-1 dPe ^ 

3— + {P2 + (^2)Pe + {P6 

p dp 



^&)P2 + [pa+UJ4: 



fin 



3 539 



+ 



17 



P2P4-P2P2P4 + ( ^- Y^sl 



19 



49 



^53 



P2PI (3.41) 



Qo] UJ2P0 + -P9P2 + -P2 \ (79 

64 9 ^ 6^^ V 64 9 ^2 



fp6,D_ 
dp^ 



P 



^"^Pe p , 1 p2 

" + 2^^^*^ + 4^^ 



1 1 p4 

4 2^ 8^ ^ 16 



1 



-P2P2PA 

11 2 3 
y 52 + 



0. 



(3.42) 



We remind the reader that the only non- vanishing lo^ for 2<Z)<6isa;2 = — 1, and we 
will show in Subsection |3.5| , that in general, for D = Q the only nonzero component is 1^4. 
The above expressions, especially those for <I> and simplify considerably for symmetric 
[/($) potentials, in which case 172 = 0. 



3.4 Free scalar field in 2 < D < 6 dimensions 

If <I> is a free massive field with potential U (<I>) = m^<I>^/2, after scaling out m and k no 
parameters remain in the equations determining Pi and pi. The spatially localized nodeless 
positive solution of the ordinary differential equations ( 3.15| ), ( 3.16| ), and the corresponding 
solution of (|]3|), (^) , ( CT) and (|1|) can be calculated numerically. For D = ?> the 
obtained curves are shown on Figs. || and |2|. The obtained central values of Pi and pi for 
i = 2, 4, 6 in D = 3, 4, 5 are collected in Table |l|. 
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Figure 1: The first three Pk functions for the free scalar field case in D = 3 spatial dimensions. 
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Figure 2: The pk functions for the free scalar field case in D = 3 dimensions. 



The chosen central values make all functions 
Pi and Pi, and consequently $j and ipi, tend to 
zero for p — oo. Although for i > 4 Pj and pi 
are not monotonically decreasing functions, their 
central values represent well the magnitude of 
these functions. Generally, the validity domain 
of an asymptotic scries ends where a higher or- 
der term starts giving larger contributions than 
previous order terms. For D = 3 the sixth order 
£ expansion can be expected to be valid even for 
as large parameter values as e = 1. For = 4 
this domain is e < 0.7, while for D = 5 it de- 
creases to e < 0.22. 





D = 3 


D = 4 


D = 5 


P2c 


2.04299 


7.08429 


28.0399 


P2c 


1.93832 


4.42976 


14.90729 


Pic 


0.658158 


-5.93174 


-348.868 


Pic 


0.686532 


-4.08270 


-200.353 


Pec 


0.557141 


27.3950 


9532.72 


P6c 


0.541339 


17.9090 


5500.18 



Table 1: Central values of the first three 
functions Pi and pi for the free scalar field 
in 3, 4 and 5 spatial dimensions. 
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3.5 order for D = 6 



P4 = ^{2Z + Z), p4+0J4=^iZ-z) 



equations (|3.36D and (|3.37| ) decouple, 



As we have already stated in Subsection 3.1, if D = 6 then u>2 = 0, s = S and the explicit 
form of the solution of the SN equations is given by (3.22). Introducing the new variables 
z and Z by 

(3.43) 

(3.44) 
(3.45) 

(3.46) 
(3.47) 



(fz 5 dz 
dp^ p dp 



sz + -Pis'^ 



d^Z 5 dZ 
dp^ p dp 



0. 



where the constants j3i and (32 are defined by the coefficients in the potential as 

80 

/3i = 1 + —gl - 8(73 , 



. 80 2 8 



For a free scalar field with U{^) = <I>^/2 we have f3i = P2 = 1- The general regular solution 
of ( f3.44 ) can be written in terms of the (complex indexed) associated Legendre function P 
as 

144/3ia^(14 + 6a2p^ + "V^) , ^"1 p2 ^1 



+ 



00 IS z„ 



13(1 + aV)^ ' aV 

where Ci is some constant. The limiting value at p 
-297.495 Ci. The regular solution of ( p^ ) is 

^'^[1 - a^p^) 



1 



(3.48) 



-Cicosh(7rV23/2)/7r 



Z 



'2OL 



7(1 + a2p2)3 
324/32aV^(220 + lOOa^pS 



ln(l + aV^) 

16aV^ 



(3.49) 



The limiting value at p 



35(1 

00 is Zoo = 



a'^p'^Y 



—324/320^/35, independently of C2. Since P4 must 
648/?2a'^/35, fixing the constant Ci. Since the mass 



tend to zero, according to (|3.43D , z^ 
of the field ^ is intended to remain m = 1, the limit of p4 also has to vanish, giving 



0)4 



324 4 
35 



(3.50) 



This expression is not enough to fix lo^ yet, since a is a free parameter. If /?2 > then it is 
reasonable to use ( ^.50 ) to set UJ4 = —1, thereby fixing the free parameter a in the order 
component of ^ and if. The change of the so far undetermined constant C2 corresponds 
to a small rescaling of the parameter a in the expression ( p.22| ). Its concrete value will 
fix the coefficient cjq in the expansion of the frequency. The homogeneous parts of the 
differential equations at higher e order will have the same structure as those for P4 and p4. 
Choosing the appropriate homogeneous solutions all higher components can be set to 
zero, yielding 

' for D = 6 if /32 > . 



UJ 



(3.51) 
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This expression is valid for the free scalar field case with potential C/($) = jl in = 6, 
since then /32 = 1- For certain potentials /32 < 0, and one can use ( |3.50D to set UJ4 = I. 
This case is quite unusual in the sense that the frequency of the oscillon state is above the 
fundamental frequency m = 1. In the very special case, when ^2 = the frequency differs 
from the fundamental frequency only in or possibly higher order terms. 

3.6 Total energy and dilaton charge of oscillons 

Substituting (|]3|) and ( ^!40|) into the expression ( p.3[) of the total energy, we get 



where 
Eo 



IT 



D/2 



vr 



D/2 



dpp^-lP2(2P4-P2) 



Since P2 = 25, for 2 < L» < 6 we can use (|334| ) and (|335| ) to get 



(3.52) 



(3.53) 



T{D/2) 



{D-2)si 



(3.54) 



The numerical values of si, Eq and Ei for D = 3, 4, 5 are listed in Table 

To the calculated order, i.e. up to e^-^, for 
D = 3 and D = 4 the energy is a monotonically 
increasing function of e, while for D = 5 there is 
an energy minimum at e = 0.2288. This result can 
only be taken as an estimate, as the validity domain 
of an asymptotic series ends when two subsequent Table 2: The numerical values of si 

Eq and £^1 in 3, 4 and 5 spatial dimen- 
sions. 





D = 3 


D = 4 


D = 5 


Sl 


3.50533 


7.69489 


10.4038 


Eo 


88.0985 


607.565 


1642.91 


El 


123.576 


2522.10 


31374.2 



terms are approximately equal. 

For D = 6 the leading order term in the total 
energy is 

E 



1927r^ 



(3.55) 



As we have already noted, for D > 6 there are no finite energy solutions. 

The leading order e dependence of the dilaton charge for 2 < D < 6 is given by 



Q = SiE 



A-D 



(3.56) 



where we used the definition ( |2.7| ), ( |3.24 ) and the relation p = er. The dilaton charge for 
the D = 6 oscillon is infinite. In higher orders in e the proportionality between the dilaton 
charge and energy is violated. 



4. Time evolution of oscillons 

In this section we employ a numerical time evolution code in order to simulate the actual 
behaviour of oscillons in the scalar-dilaton theory. We use a fourth order method of line 
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code with spatial compactification in order to investigate spherically symmetric fields |46|. 
Our aim is to find configurations which are as closely periodic as possible. To achieve this, 
we use initial data obtained from the leading terms of the small amplitude expansion 
(3.39) and ( |3.40| ). The smaller the chosen e is, the more closely periodic the resulting 
oscillating state becomes. However, for moderate values of e, it is possible to improve the 
initial data by simply multiplying it by some overall factor very close to 1. 

The main characteristics of the evolution of small amplitude initial data depend on 
the number of spatial dimensions D. For D = 3 and D = 4 oscillons appear to be stable. 
If there is some moderate error in the initial data, it will still evolve into an extremely long 
living oscillating configuration, but its amplitude and frequency will oscillate with a low 
frequency modulation. We employ a fine-tuning procedure to minimize this modulation 
by multiplying the initial data with some empirical factor. For D = 5 and D = 6 small 
amplitude oscillons are not stable, having a single decay mode. In this case we can use the 
fine-tuning method to suppress this decay mode, and make long living oscillon states with 
well defined amplitude and frequency. Without tuning in D = 5 and D = 6, in general, 
an initial data evolves into a decaying state. The tuning becomes possible because there 
are two possible ways of decay. One with a steady outwards flux of energy, the other is 
through collapsing to a central region first. 

Having calculated several closely periodic oscillon configurations, it is instructive to 
see how closely their total energy follow the expressions (|3.52| )-( 3.55 ). Apart from checking 
the consistency of the small amplitude and the time-evolution approaches, this also gives 
information on how large e values the small amplitude expansion remains valid. The 
parameter e for the evolving oscillon is calculated from the numerically measured frequency 
by the expression e = \/l — uj'^. The results for D = 3 are presented on Fig. |3|. In contrast 
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Figure 3: Total energy of three-dimensional oscillons as a function of the parameter e. 

to general relativistic oscillatons, there is no maximum on the energy curve. This indicates 
that all three dimensional oscillons in the dilaton theory are stable. 

The e dependence of the energy for D = 5 is presented on Fig. ^ There is an energy 
minimum of the numerically obtained states, approximately at e = 0.21, above which 
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Figure 4: Total energy of five-dimensional oscillons as a function of the parameter e. The vertical 
line at e = 0.21 shows the place of the energy minimmn. States to the right of it are stable, while 
those to the left have a single decay mode. 



oscillons are stable. The place of the minimum agrees quite well with the value e = 0.2288 
calculated in Subsection [s]^ using the first two terms of the small amplitude expansion. 
The behaviour of the energy close to the minimum is shown on Fig. ^. 
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Figure 5: The region of Fig. ||near the energy minimum. 

We have also constructed oscillon states for D = 6 dimensions. These oscillons have 
quite large energy, due to the slow spatial decay of the functions ^ and ip. For free massive 
scalar fields, oscillons have frequency given by ( ^.51 ), i.e. an initial data with a given e 
value will evolve to an oscillon state with frequency approximately following u> = \/l — e'^. 
However, there are potentials, for which the oscillation frequency is above the threshold 
Lu = 1. For example, this happens for the potential C/($) = ^^{^ — 2)^/8 with the choice 
K = 1/2. 

In conclusion, in the dilaton-scalar theory oscillons follow the stability pattern observed 
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in the case of self-interacting scalar and Einstein-Klein-Gordon theory; if e decreases with 
decreasing energy, oscillons are stable, while if e increases with decreasing energy, oscillons 
are unstable. In other words if the time evolution (i.e. energy loss) of an oscillon leads 
to spreading of the core, the oscillon is stable, while oscillons are unstable, if they have to 
contract with time evolution. The decreasing or increasing nature of the energy, and hence 
empirically the stability of the oscillating configurations, is well described by the first two 
terms of the small amplitude expansion ( 3.52| ). The result following from Eq. ( 3.52| ) shows 
the existence of an energy minimum for D > 4. This provides an analytical argument for 
the existence of at least one unstable mode. In particular, for D = 5 spatial dimensions the 
frequency separating the stable and unstable domains is determined by the small amplitude 
expansion to satisfactory precision. 

In order to study the instability in more detail numerically, we compared the evolution 
of two almost identical initial data obtained from the small amplitude expansion with 
e = 0.05. In order to make the unstable state long living, a fine tuning procedure is 
applied, multiplying the amplitude of the initial data by a factor with value close to 1.0178. 
The multiplicative factors used in the two chosen initial values differ by 2.2 x 10^^^. One 
of the two initial data develops into a configuration decaying with a uniform outward 
current of energy, the other through collapsing to a high density state first. On Fig. ^ 
the time evolution of the difference of the central value of the dilaton fields in the two 
states Aip = ^pi — (p2 is shown. The curve follows extremely well the exponential increase 




Figure 6: Increase in the difference of the dilaton field (p for two similar configurations. 



described by 

Aif = 6.583 X 10"^^ exp(0.003157t) , (4.1) 

showing that there is a single decay mode growing exponentially. The difference of the 
scalar fields, A<I> = <I>i — ^2, grows with the same exponent. The spatial dependence of 
the decaying mode is illustrated on Fig. 0, where A$ is plotted at several moments of time 
corresponding to the maximum of <I>i at the center. 
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Figure 7: Radial behaviour of the difference of the scalar fields of two very similar configurations. 
At the chosen moments of time the scalar is maximal at the center, and the subsequent moments 
are separated by 32 oscillations. 



Our numerical results strongly indicate that oscillons in the scalar-dilaton theory are 
unstable for D > 4, and they admit a single decay mode. For the single scalar field 
system the instability arises for D > 2 (see Ref. [31 1), but the decay modes have been 
calculated analytically only in some very special cases. The scalar theory with potential 
(/{(p) = cf^ll — Incf^) admits exactly time-periodic breathers in any dimensions. The 
stability of these breathers in three dimensions has been investigated in detail in [47||48|. 
It has been found that these breathers always admit a single unstable mode. It needs 
further studies whether an analysis along the lines of Ref. [^] can also be applied to more 
general potentials in the small amplitude limit, and whether it can be generalized to the 
case when the scalar is coupled to a dilaton field. 



5. Determination of the energy loss rate 

Although oscillons are extremely long living, generally they are not exactly periodic. In 
this section we calculate how the energy loss rate depends on the oscillon frequency for 
small amplitude configurations. To simplify the expressions in this section we consider a 
massive free scalar field, i.e. C/($) = m^$^/2. We assume that 2 < D < 6, since then the 
scalar field tends to zero exponentially for large p. 

The outgoing radiation will dominantly be in the dilaton field and the radiation am- 
plitude will have the e dependence: e ex.p {—2Q d / e) . In Refs. [30| and |31| we have used 



two different methods for determining the e independent part of the radiation amplitude: 
Borel summation and solution of the complexified mode equations numerically. In this 
paper we will use an analytic method based on Borel-summing the asymptotic series in the 
neighborhood of its singularity in the complex plane. Other potentials which are symmet- 
ric around their minima can be treated analogously. If the potential is asymmetric only 
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the numerical method could be used. This phenomenon is in complete analogy with the 
problem arising with a single scalar field considered in Ref. pO[| . 



5.1 Singularity of the small e expansion 



We first investigate the complex extension of the functions obtained by the small amplitude 
expansion in Sec. ^. Extending the solutions s and S of the Schrodinger-Newton equations 
(3.19) and (3.20) to complex p coordinates they both have pole singularities on the imag- 



inary axis of the complex plane. We consider the closest pair of singularities to the real 
axis, since these will give the dominant contribution to the energy loss. They are located 
at p = ±iQD- The numerically calculated values oi Qd are listed in Table ^ for spatial 
dimensions D = 3, 4, 5. 

Let us measure distances from the upper singularity by a 
coordinate R defined as 



P = iQo + R- 



(5.1) 



Close to the pole we can expand the SN equations, and obtain 
that s and S have essentially the same behaviour. 



D 


Qd 


3 


3.97736 


4 


2.30468 


5 


1.23595 



±s 



_6_ 
R^ 



6i{D-l) {D-l){D-5l) 



5QdR 



50Q2 



+ OiR), (5.2) 



D 



Table 3: The distance 
Qd between the real axis 
and the pole of the funda- 
mental solution of the SN 
equation for various spa- 
even though they clearly differ on the real axis. Since for sym- ^ial dimensions D. 
metric potentials we can always substitute ^> by — <I>, we choose the positive sign for S 
in (5.2). This choice is compatible with the sign of S used on the real axis at the small 
amplitude expansion section. We note that for D > 1 there are logarithmic terms in the 



expansion of s and S, starting with terms proportional to R'^lnR. According to ( 3.1^ 
the functions determining the leading parts of (p and <I> in this case are 



P2 = s + 1, 



P2 = 2S. 



(5.3) 



Substituting these into the equations ( ^.36 ) and ( ^.37|) , the order contributions p4 and 
P4 can also be expanded around the pole 



P4 

P4 - 2p4 



1161 524i{D-l)lnR c_3 .^(^^^R 
~52^+ 35QdR^ 

81 18i{D - 1) / 1 
13R^ 5QdR^ \R^ 



(5.4) 
(5.5) 



where the constant c_3 can only be determined from the specific behaviour of the functions 
on the real axis, namely from the requirement of the exponential decay of P4 for large real 
P- 

5.2 Fourier mode expansion 



Since all terms of the small amplitude expansion (3T) are asymptotically decaying, i.e. 
localized functions, the small amplitude expansion can be successfully applied to the core 
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region of oscillons. However it cannot describe the exponentially small radiative tail re- 
sponsible for the energy loss. Instead of studying a slowly varying frequency radiating 
oscillon configuration it is simpler to consider exactly periodic solutions having a large 
core and a very small amplitude standing wave tail. We look for periodic solutions with 
frequency uj by Fourier expanding the scalar and dilaton field as 



Nf 



$ = ^ ^Iffc cos{ku}t) , 



k=0 



Nf 

(f = Vfc cos{kujt) . 

k=0 



(5.6) 



Although, in principle, the Fourier truncation order Np should tend to infinity, one can 
expect very good approximation for moderate values of Np. In (^) we denoted the 
Fourier components by psi instead of phi to distinguish them from the small e expansion 



components in (3.1). Since in this section we only deal with an self-interaction free scalar 
field with a trivially symmetric potential. 



2k 



0, 



■02^+1 = , for integer k . 



(5.7) 



We note that the absence of sine terms in ( |5.6D is equivalent to the assumption of time 
reflexion symmetry at t = 0. This assumption appears reasonable physically, and we have 
seen in Sec. ^ that it holds in the small amplitude expansion framework. 

For small amplitude configurations we can establish the connection between the ex- 
pansions (|3.lD and (|5.6|) by comparing to (|3.39| ) and ( |3.4C ), obtaining 



^-1 

'4'2 



256 ^ ^ 



p2 
.4 -'^2 

' 16 



32 



- iP2 - l)Pi 



dR 



(5.8) 
(5.9) 
(5.10) 

(5.11) 



Let us define a coordinate y for an "inner region" by R = ey. This coordinate will 
have the same scale as the original radial coordinate r, since they are related as 



+ y- 



(5.12) 



The "inner region" \R\ ^ 1 is not small in the y coordinate; if e ^ then e\y\ = \R\ <ti 1 
but \y\ — > 00. Using the coordinate y and substituting (^)-( [5l^ ) into the small amplitude 
Fourier mode expressions ( |5.8D -( 5Tl| ), we obtain that the leading asymptotic behaviour of 
the Fourier modes for |y| — > cxd can be written as 



^-1 



12 
I2 



999 



+e 



26y4 
6i{D - 1 



+ elne 



6m{D - 1) 



27 81i(Z) 
%6 



— + 



35QDy^ 
108 In y 



7y2 



2 + 



2c_ 



+ 



1) 



20Qz5y5 



+ ... 



(5.13) 
(5.14) 
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1161 , 324im 



1) 



6i(L' - 1) /541ny 



9 

V'2 = 4 

y 



^Qdv 
m{D 



7y2 



1 + 



C-3 



+ .. 



11 



These expressions are simultaneous series in 1/y and in e. 



(5.15) 
(5.16) 



5.3 Fourier mode equations 

In order to obtain finite number of Fourier mode equations witli finite number of terms, 



wlien substituting (5^) into tlie field equations (|2.5|) and (|2.6D we Taylor expand and 
truncate the exponential 



(5.17) 



fc=0 



We need to carefully check how large iVg should be chosen to have only a negligible influence 
to the calculated results. For n < Np the Fourier mode equations have the form 

d2 



dr^ 



+ 



D-1 d 
dr 



+ n^uj'^ - 1 1 ^'r 



dr^ 



D-1 d 2 2,, 

H — + n uj ] tpn 

r dr 



fn 



(5.18) 
(5.19) 



where we have collected the nonlinear terms to the right hand sides, and denoted them with 
Fn and fn- These are polynomial expressions involving various and tpk, with quickly 
increasing complexity when increasing the truncation orders Np and N^. The solution of 
(|5.18 ) and ( 5.191) yields the intended quasibreathers, with a localized core and a very small 
amplitude oscillating tail. For small amplitude configurations the functions and ipk will 
have poles at the complex r plane, just as we have seen in the small amplitude expansion 
formalism. In order to calculate the tail amplitude it is necessary to investigate the Fourier 
mode equations instead of the equations obtained in Sec. ^. Although in the Fourier 
decomposition method we have not defined a small amplitude parameter yet, motivated 
by ( |3.30| ), we can, in general, define e as 



(5.20) 



Dropping 0{e'^) terms, in the neighborhood of the singularity the mode equations take the 
form 



d^ D-1 d 

dy2 iQ^ dy 



f d D-l d 2 , , 

\dy^ iQd dy 



(5.21) 
(5.22) 



We look for solutions of these equations that satisfy ( 5.13| )- (|5.16 ) as boundary conditions 
for \y\ — > 00 for — 7r/2 < argy < 0. This corresponds to the requirement that the functions 
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decay to zero without any oscillating tails for large r on the real axis. The small correction 
corresponding to the nonperturbative tail of the quasibreather will arise in the imaginary 
part of the functions on the Rey = axis. 

5.4 e ^ limit near the pole 

For very small e values one can neglect the terms proportional e on the left hand sides 
of ( 5.21 ) and ( ^.22|) . In this limit the there is no dependence on the number of spatial 



dimensions D. We investigate this simpler system first, and consider finite but small e 
corrections later as perturbations to it. We expand the solution of ( ^.21| ) and ( |5.22| ) (with 
e = 0) in even powers of 1/y, 



CO ^ oo ^ 

j=k+l ^ j=k+l ^ 

We illustrate our method by a minimal system where radiation loss can be studied, namely 
the case with Np = 3 and = 1. Then the mode equations are still short enough to 
print: 

^ = ^{^Po - + ^1) + ^V'2^'i(^i + 2^3) , (5.24) 

__i = _ _^2(^i + ^-3) , (5.25) 

^ + 4^2 = 7^'i(V'o - l)(^i + 2^-3) + IM'^I + ^1^3 + ^3) > (5-26) 
dy^ 4 4 

^ + 8^3 = -2^2^-1 -V'o^s. (5.27) 



When looking for solution of these equations in the form of the 1/y^ expansion ( [5.23| ) , only 
one ambiguity arises, the sigr 
the expansion turn out to be 



one ambiguity arises, the sign of A'^\ Choosing it to be negative, the first few terms of 



i'0 = ---^,+O{^] , (5.28) 
*i = -^-£^ + 0(^) > (5-29) 



6 


837 


y2 


52y4 ^ 


12 


459 


y2 


26y4 ^ 


9 


1845 


¥~ 


52y6 ^ 


27 


2565 


4^ " 


416y^ 



^2 = -i^-^ + 0[^j , (5.30) 
*3 = -^-^ + of^) . (5.31) 



The first terms agree with those of ( 5.13 )-( 5.16 ) obtained by the small amplitude expansion. 
The difference in the 1/y^ terms of ipo and ^1 are caused by the too low truncation for the 
Taylor expansion of the exponential. For >2 these terms agree as well. 

When increasing Np and growing number of additional terms appear on the right 



hand sides of ( 5.24D -( p.27 ), and the number of mode equations rise to A'^;' + 1. These 
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complicated mode equations can be calculated and expanded using an algebraic ma- 
nipulation program. However, apart from a factor, the leading order behaviour of the 



coefficients a^^ and AX'"' for large n will remain the same as that of the minimal system 



(5.24)-( p.27 ). The large n behaviour of these coefficients will be essential for the calculation 
of the nonperturbative effects resulting in radiation loss for oscillons. 

Starting from the free system, consisting of the linear terms on the left hand sides, it 
is easy to see that the mode equations are consistent with the following asymptotic (large 
n) behavior of the coefficients, 



^0 ' ^1 



a. 



in) 
2 

(n) 



(2n - 1) 
4" 



« 4") 



(5.32) 
(5.33) 



where k is some constant. The value of k can be obtained to a satisfactory precision 
by substituting the 1/y expansion into the mode equations and explicitly calculating the 
coefficients to up to high orders in n. In practice, using an algebraic manipulation software, 
we have calculated coefficients up to order n = 50. The dependence of k on the order of 
the Fourier expansion is given in Table 

The results strongly indicate that in the Nj ,Ne^oo 
limit k = 0. We do not yet understand what is the deeper 
reason or symmetry behind this. Hence, instead of ( 5.32| ) 
and ( ^.33| ), the correct asymptotic behavior is 



in) 



a, 



in) j^in) ^(n) 







On 



Ki-1) 

« 4") . 



(2n- 1)! 



(5.34) 
(5.35) 







k 


3 


5 


3.71 X 10" 


-3 


4 


7 


3.12 X 10" 


-6 


5 


8 


6.03 X 10- 


-9 


6 


10 


4.61 X 10" 


13 



Table 4: Dependence of the con- 
stant k on the considered Fourier 
components Np- The second col- 
umn lists the minimal exponential 
expansion order which is nec- 
essary to get the k value with the 
given precision. 



Taking at least Np = 6 and N^, = 9, the numerical value 
of the constant turns out to he K = —0.57 it 0.01. The 
above results indicate that the outgoing radiation is in 
the ^'3 scalar mode instead of being in the ip2 dilaton 
mode. This conclusion is valid only in the framework of 
the approximation employed in the present subsection, 

i.e. when dropping the terms proportional to e in (|]2lD and ( |532D . As we will see in the 
next subsection, the situation will change to be just the opposite when taking into account 
£ corrections. 

All terms of the expansion ( 5.23| ) are real on the imaginary axis Rey = 0. However, 
using the Borel-summation procedure it is possible to calculate there an exponentially 
small correction to the imaginary part. We will only sketch how the summation is done, 
for details see |30] and |37|. We illustrate the method by applying it to ^'3. The first step 
is to define a Borel summed series by 



Viz) 



^ A 



in) 



n=2 



(2n)! 



00 
n=2 



-ir 

2n V\/8y 



2n 



2 



1 + ^ 



(5.36) 
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This series has logarithmic singularities at z = ibi-v/8. The Laplace transform of V{z) will 
give us the Borel summed series of ^'3(2/) which we denote by ^3(1/) 



00 /J. 



= dte-'Vi^-j . (5.37) 

The choice of integration contour corresponds to the requirement of exponential decay on 
the real axis. The logarithmic singularity of V {t/y) does not contribute to the integral 
and integrating on the branch cut starting from it yields the imaginary part 

Irn^-M = dte-'— = — exp (-iVSy) . (5.38) 

JiVSy 2 2 \ J 

A similar calculation for the "02 dilaton mode yields 

Im^2(y) = exp (-2%) . (5.39) 

Since A: = 0, this mode is vanishing now. However, as we will show in the next subsection, 
when taking into account order e corrections a similar expression for 1(^2 with exp {—2iy) be- 
haviour arise, which, due to its slower decay, will become dominant when Imy — > —00. The 
continuation to the real axis of these imaginary corrections turns out to be closely related 
to the asymptotically oscillating mode responsible for the slow energy loss of oscillons. 

5.5 Order e corrections near the pole 

Before discussing the issue of matching the imaginary correction calculated in the neigh- 
borhood of the singularity to the solution of the field equation on the real axis we deal with 
the corrections arising when taking into account the terms proportional to e in the mode 



equations ( ^.21 ) and ( |5.22| ). We denote the solutions obtained in the previous subsection 



by il^n^ and '^n \ and linearize the mode equations around them by defining 

^Pn = V'i"^ +^n, ^-n = + • (5.40) 

The mode equations take the form 

/ 2\r ^ D-l d4°' ^»fn7 J, K .,1 



/ 2 .At D-l d^-r ^ dFn~ V- dFn ~ 



where the partial derivatives on the right hand sides are taken at = '^n^ and ipn = V'r?^ • 
The small dimensional corrections V'n and ^„ have parts of order both elne and e. 

The linearized equations ( 5.41| ) and ( 5.42| ) are solved to e\ne order by the following 
functions: 

d^^(°^ 

V'„ = elneC^^, (5.43) 
ay 

d^l°) 

^'„ = elneC— ^, (5.44) 
Ay 
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where C is an arbitrary constant. The reason for this is quite simple: in elne order the 
terms proportional to e on the left hand sides are negligible and we get the e = equation 
linearized about the original solution. Our formula simply gives the zero mode of this 
equation. The constant C is determined by the appropriate behaviour when continuing 
back our functions to the real axis. This can be ensured by requiring agreement with the 
first few terms of the small amplitude expansion formulae ( 5.13| )-( p.l6| ), yielding 



C 



27i{D - 1) 



35Q 



D 



(5.45) 



In the small amplitude expansion ( |5.13| )-( |5.16| ) to every term of order e In e corresponds 
a term of order e which we get by changing Ine to Iny. Thus, we define the new variables 
^l^J^ and to describe the e order small perturbations by 



V'n = £ In £ C " + e C In y " + tp^^ 
dy \ dy 



d^ 



(0) 



= elneC—^ +e\ Clny , 
dy \ dy 



(5.46) 
(5.47) 



Substituting into the linearized equations ( 5.41 ) and ( ^.42 ) we see that all terms containing 
Iny cancel out, 



dy2 



C 



+ n'j Vn + 
dfn -J 



2y 



dy2 



(0)^ 



dy 



+ 



D-ldip, 



(0) 



fd' 2 



C_ 



2y 



d^ 
dy2 



d^ 

dy 



iQd dy 



+ 



D 



1 d^f 



iQd dy 



(5.48) 



(5.49) 



If C is given by ( ^.45 ), -i/'n turn out to be algebraic asymptotic series which are 

analytic in y. Let us write their expansion explicitly: 



2k+l 



-°2fc+l„,2n-l 



2k 



E 



n=k+l 



Substituting these and the expansions (|5.23|) for ipn and into ( ^.481) and ( 5.49 ), it is 



n=k+l 

(0) 



2k y2n-l 



(5.50) 



possible to solve for the coefficients hl"^ and BV'' , up to one free parameter. Comparing 



(n) 



to (|5.13D-(B.16D it is natural to choose this free parameter to be 6n = c_3. Similarly to 
^ CI) 

that case, 6q will only be determined by the requirement that the extension to the real 
axis represent a localized solution. Furthermore, leaving C a free constant and requiring 
the absence of logarithmic terms in the expansion of V'fc ^-nd ^'^ yields exactly the value of 
C given in ( 5.45| ). 
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Eq. ( |5.48| ) is consistent with the asymptotics 



62 ~ ikd (-1) —TT, 



2)\ 



22n-l 



(5.51) 



where kfj is some constant. Since the leading order result for is given by ( 5.34 ), if 
kf) 7^ 0, the co( 
the value of k£ 
1/y, obtaining 



7^ 0, the coefficients follow the hierarchy 62"^ ^ ^3" ■ In order to be able to extract 
the value of kf) we have calculated by solving the mode equations to high orders in 



kD = 1M0^^—^ . (5.52) 
Qd 

The displayed four digits precision for ki) can be relatively easily obtained by setting 

(n) 

Np > 4:, Ne > 5 and calculating bX to orders n > 25. We note that there is also a term 

(2) (n) 

proportional to the unknown 6g = c_3 in each 62 > giving a c_3 dependent ho- Luckily, 
the influence of this term to k^i quickly becomes negligible as Np and A'^e grow, making 
the concrete value of c_3 irrelevant for our purpose. 

The Borel summation procedure can be done similarly as in Eqs. (|5^-(p8|). On the 



imaginary axis 1^2 is real to every order in 1/y, however it gets a small imaginary correction 
from the summation procedure given by 

ImV^2(y) = e exp (-2%) . (5.53) 
5.6 Extension to the real axis 



Solutions of the Fourier mode equations ( 5.1S| ) and ( 5.1S| ) can be considered to be the sum 



of two parts. The first part corresponds to the result of the small amplitude expansion, 
the second to an exponentially small correction to it. The small amplitude expansion 
is an asymptotic expansion, it gives better and better approximation until reaching an 
optimal order, but higher terms give increasingly divergent results. The smaller e is, the 
higher the optimal truncation order becomes, and the precision also improves. The small 
amplitude expansion procedure gives time-periodic localized regular functions to all orders, 
characterizing the core part of the quasibreather. Their extension to the complex plane is 
real on the imaginary axis. Furthermore, the functions obtained by the e expansion are 
smooth on large scales, missing an oscillating tail and short wavelength oscillations in the 
core region. On the imaginary axis, to a very good approximation, the small second part 
of the solution of the mode equations ( 5.18| ) and ( 5.19| ) is pure imaginary, and satisfies 



the homogeneous linear equations obtained by keeping only the left hand sides of these 
equations, because the quasibreather is a small-amplitude one. In the "inner region" it is 
of order , while on the real axis its amplitude is of order , hence to leading order the 
quasibreather core background does not contribute. In the previous subsection we have 
determined the behaviour of this small correction close to the poles. Now we extend it to 
the real axis. 

In the "inner region", close to the pole, the function Im^2 given by ( ^.531) solves 
the homogeneous linear differential equations given by the left hand side of (|5^ ). The 
extension of this function to the real axis will provide the small correction to the small 
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amplitude result mentioned in the previous paragraph. We intend to find the solution 1^2 
of the left hand side of ( 5.19 ), which reduces to the value given by ( ^.53| ) close to the upper 
pole, where r = iQz)/e + y, and behaves as 



Im ■02 (y) 



-e ^ exp (2iy) . 



(5.54) 



near the lower pole, where r = —iQc/s + y- We follow the procedure detailed in |31]. The 
resulting function for large r is 



^^2 



ikoTT ( Qd 



er 



exp 



exp(-2ir) - (-i)(^-i)/2 exp(2ir) 



(5.55) 

The general solution of the left hand side of ( 5.19| ) can be written as a sum involving Bessel 
functions Jn and 1^, which have the asymptotic behaviour 



/ 2 / VTT ■K\ 

J.(x)-y-cos(:r----j, 



for X +00. The solution satisfying the asymptotics given by ( 5.55| ) is 



^2 = V^J^Yr,/2-i{2r), 



where the amplitude at large r is given by 



(D-l)/2 



exp 



2Q 



D 



(5.56) 
(5.57) 

(5.58) 
(5.59) 



For D > 2 the function given by ( 5.5§| ) is singular at the center, due to the usual central 
singularity of spherical waves. Since the amplitude of the quasibreather core is proportional 
to e^, and its size to 1/e, for small e it is possible to extend the function -02 in its form ( 5.5^ ) 
to the real axis into a region which is outside the domain where ■02 gets large, but which 
is still close to the center when considering the enlarged size of the quasibreather core. 
When extending this function further out along the real r axis, because of the large size of 
the quasibreather core, the nonlinear source terms on the right hand side of ( 5.19| ) are not 
negligible anymore, and the expression ( 5.5§| ) for -02 cannot be used. What actually happens 
is that ^02 tends to zero exponentially as r ^ cxd. This follows from the special choice of 
the "inner solution" close to the singularity; namely, we were looking for a solution which 
agreed with the small amplitude expansion for Rey — > 00. The small amplitude expansion 
gives exponentially localized functions to each order and we also required decay beyond all 
orders when choosing the contour of integration in the Borel summation procedure. 

By the above procedure we have constructed a solution of the mode equations which is 
singular at r = 0. The singularity is the consequence of the initial assumption of exponential 
decay for large r. The asymptotic decay induces an oscillation given by (|]5D in the 
intermediate core, and a singularity at the center. In contrast, the quasibreather solution 
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has a regular center, but contains a minimal amplitude standing wave tail asymptotically. 
Considering the left hand side of ( p.l9| ) as an equation describing perturbation around 
the asymptotically decaying solution, we just have to add a solution 6ip2 determined by 
the amplitude ( ^.59| ) with the opposite sign of ( |5.58| ) to cancel the oscillation and the 
singularity in the core. This way one obtains the regular quasibreather solution, whose 
minimal amplitude standing wave tail is given as 



-V^Z^^D/2-i{2r)cos{2t) (5.60) 



r{D-i)/2 



sm 



2r-{D- 1)^ cos(2t). (5.61) 



Adding the regular solution, where Y is replaced by J, would necessarily increase the 
asymptotic amplitude. 

If we subtract the incoming radiation from a QB and cut the remaining tail at large 
distances, we obtain an oscillon state to a good approximation. Subtracting the regular 
solution with a phase shift in time, we cancel the incoming radiating component, and obtain 
the radiative tail of the oscillon, 

(pose = -^-^HZl [^^D/2-i(2r)cos(2t) - J^/2_i(2r) sin(2t)] (5.62) 



a 



D 



(D-l)/2 



sm 



r 



2r - (Z) - 1)^ - 2t . (5.63) 



The radiation law of the oscillon is easily obtained now, 

where the constant ko is given by ( |5.52| ). If we assume adiabatic time evolution of the 
e parameter determining the oscillon state, using Eqs. (|]5|) and giving £^ as a 



function of e, we get a closed evolution equation for small amplitude oscillons, determining 
their energy as the function of time. 

For the physically most interesting case, = 3 we write the evolution equation for e 
and its leading order late time behavior explicitly: 

de ( 15.909 \ 

— = -30.29 exp [ — j (5.65) 

15.909 ^ 1401.6 , _ 

£ « — , E — . (5.66) 
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